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ABSTRACT 


Measurements of the ratio of the elastic electron- 
deuteron scattering cross section to the elastic electron- 
proton scattering cross section were made for low momentum 


transfers. Augmented by the work of Topping [5], these data 


2 2 2 


Span a range of q~ from 0.05 fm “ to 0.50 fm “~. From these, 


the electric form factor of the neutron, Gr , was extracted 
n 


using three Lomon-Feshbach wave functions for the deuteron, 


Z 
ne versus q 
n 


were measured to be 0.0218 + 0.0022 for LFl, 0.0220 + 0.0022 


(LF1, LF5, and LF15). Values of the slope of G 


feet >, and 0.0185 + 0.0022 for LF15. Although the data 
seems to prefer model LF15, all were found to be in agree~ 
ment with the thermal neutron-electron slope of 0.0193 + 
0.0004. The root mean square structure radius of the 
deuteron was also extracted, using a method devised by 
Schumacher, and was found to be 1.9665 + 0.0045 fm. Rela- 
tivistic corrections proposed by Casper and Gross [6, 18] 
were Peipares With the relativistic correction of Friar [8] 
and no significant differences were found in the region 


investigated. 
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iL. TNEROpUC TION 


The deuteron being composed of one proton and one 
neutron is the only stable two body nucleus and, as such, 1s 
the simplest of the nuclei except for the hydrogen nucleus. 
The proton exists as a stable particle in the form of a 
hydrogen nucleus, but the isolated neutron is unstable. In 
order to obtain information concerning the neutron, this 
work will compare electron scattering from the deuteron and 
proton and extract the electric form factors of the neutron, 


Gp , Which gives information on the charge structure of the 
n 
neutron. 


0) 
This experiment measured G at various values of q 


FE; 
n 


(the square of the Lorentz invariant four-momentum trans- 


ferred from the scattered electron to the target nucleus) 


2 2 


mrioging from 0.4 fm ~ to 0.5 fm ~, in order to resolve 


differences among several earlier experiments [1,2,3,4]. 
The results of this work augmented with those of Topping [5] 


provide values of the neutron charge form factor for the 


2 Z 


momentum range q? =somg> fm < to O850 fm . “By comparing 


these Gp which are obtained with the use of several dif- 
n 
ferent wave functions, to the G,, predicted by the very 
n 


accurate thermal neutron experiment of Krohn and Ringo [4], 
a model that best represents the structure of the deuteron 
may be found. Three different models containing D-state 


Gentributions of 4.57, 5.20 and 7.55% were tested and 





although the results are not conclusive, the data favors 
bee model with /.55% D-state. Two different relativistic 
corrections were investigated, the first developed by 
Gross [6] and modified by Schumacher [7] and the second 
developed by Friar [8]. Results using both methods are 
reported here. For the momentum transfers of this experi- 
ment, it is shown that the differences between the two 
methods are negligible. From the slopes of Gp LO “Chie 
various deuteron models used, a mean Were eer sae of 
the deuteron is extracted employing a method developed by 
Schumacher [7]. 

The organization of this thesis is as follows: the 
concepts of electron scattering and the definitions of 
terms needed to describe the scattering from protons and 
deuterons are contained in Chapter II. 

Chapter III discusses briefly the experimental equipment 
and the procedures followed in order to obtain the data. 

The contents of Chapter IV discusses the methods employed 
for data peadecion. Included are discussions of corrections 
and normalizations made to the data, radiative and ioniza- 
tion corrections to the scattering spectra, and a detailed 
analysis of the sources of errors. 

Chapter V presents the results using several different 
wave functions and both methods for incorporating the 
MeeativVistic correction. Plots for G,, versus q? are 
included for all cases. Also etuGee 4 a Rosenbluth Plot 
aS an experimental verification of the conventional scaling 


law assumptions of Chapter II. 





Chapter VI presents the conclusions of this work and 
recommendations by the authors for future experiments. Five 
appendices amplify details of several topics discussed in 
this thesis. The last appendix is a compilation of the 


FORTRAN codes used in the data analysis. 


10 





Le Ory, 


A. NEUTRON CHARGE FORM FACTOR 
The cross section for an extremely relativistic electron 
scattering from a spinless point nucleus with no magnetic 


moment has been calculated by Mott [9], 








do Ze- : cos“ 5 1 
ome? — \an — 1 OF AG "ie 4 (TT=1) 
Mott 1 sin 5 
where 
2E. 
E=1+—, sie a 
Mc 
n 


is the recoil factor. 


The experimental cross section is related to the Mott 


cross section by 


(55) 2 
AN fi 
exp 


where G(q*) is the total form factor of the nucleus, which 


(S5) G* (q*) ; 
Mott 


° 


in the case of charge scatterings is a function only of the 
square of the four-momentum transfer. The form factor takes 
into account the finite extent of the nucleus. 

A theoretical expression for the experimental cross 
section for a proton has been developed by Rosenbluth [10] 
and is given in terms of the Mott cross section and total 


merm factor Ge as developed by Hand et ele Go 


11 





where oF (8) is the proton Mott cross section and 


G2 a as G G2 
a i 26 2 
G = ad i 2tT tan <=G ; 
Pp 1+ T 2 mM 
where G. and Gu are the electric and magnetic form factors 
p Pp 2 
of the proton respectively and T = a. : 
4M 
p 


mimorineci ple, Equations II-l and II-2 may be used to 


measure Gy and Gy by doing experiments at constant q with 


p Pp 
various values of incident energy and scattering angle. 


However, for low values of gq, the magnetic term is small (at 
most 15% of the total) so that it is more convenient M6 Cal 
culate the magnetic scattering using the scaling law developed 
by Lehmann et al. [12] and others [11], which was shown to 

be valid at low momentum transfers by Brooks and Sessler [13]. 
In terms of the proton monte moment, He the scaling law 


is given by 


G =u G : (II-3) 


Utilizing Equation II-3, Equation II-2 can be simplified 


to 
C2 
ae. lane 2 26 
ce ee fl + TH {1+2(1+tT) tan x} |. (II-4) 


Since the TUS term in the brackets is due only to magnetic 


Z See equation II-4 


mateco~rs, defining C_ = tus {1+2(1+t) tan 5 


p 
can be written 


12 





Gr 
ee = S82 cei 


Rewriting this and defining a proton magnetic correction 


renee , an expression for ee 
1+C,, ED 





term, Ke = is mtound; 


p = 
Ku : (IDI 3) 


From Equation II~2, it can be seen that Ge can be 


obtained from the ratio of the experimental cross section 


Pemene Mott cross section and, G2 SCanewbeouwT 1 teen in terums 





E 
Pp 
of measurable quantities. 
doP 
eo see Ky (II-6) 
p of, (8) 


An analogous method is used to develop an expression for 


Ge , the electric form factor of the deuteron. Jankus has 


d 
developed an expression for the cross section for electron 


Peatecering from the deuteron [14]. Written in a form that 
includes form factors to account for Enon eoeexrecnt Ga 


the nucleons, Jankus' results [15] are 


= _ .d Zz 2 Bo 2a 2 Dae D 2 at 

dQ ou (8) {aa (q I+ 3 Te (q“) +3 Gi, (q™){1+2(14+n)tan™ 5]7, 
ol d d 

(II-7) 

= 2 2 2 

where nq = = and G_ (q ), G. (q'), G, (q ) are respectively 
Z E O M 
4M 5 d d d 


the charge, quadrupole and magnetic moment form factors of 


the deuteron. They are given by: 


a = (G., SLE Gy ) Ch ’ G = (G pomGie) OC A 4 
Pp 7) 


13 





{(G,, a G,, ) Cy + (Gi + Gu ) Ce |. 
d p Pp n p n 
Here, Ch Pomenewechange Structure factor, C 


Bele structure factor, C 


PS the -quadru- 


Q 


L ae¢ecolmes foun the Magnetic contri- 


bution to scattering arising from charge convection in the 


deuteron and C, explains the contribution to scattering by 


5 
the magnetic moments of the proton and neutron. These four 
factors Cre Cor Che and Co are all model dependent in that 


they arise from integrals involving S and D state wave 


Penerions Of the deuteron [5]. Values for C_, C Cc and 


E Os as 
Co for momentum transfers studied were supplied by Lomon [16] 
for the Lomon-Feshbach Models LFl, LF5 and LF15 [17] which 
Sentalin percent D=-state contributions of 4.57, 5.20, and 
7.55 respectively. 


If u(r) and w(r) are the normalized S and D-state wave 


meametions, then Ch is given by 


oO 


Ce eI [u? (x) + w°(r)] 5, (4) ar , 


where Is 1S a spherical Bessel function, corresponding ex»vres- 


sions for Cor Che and Co avec counGduein Rer.  | 39). in 
the static limit (q=0), the above functions are given by 
& @) ai; Gy) (0) = MSO, 
d d 


M 
G. (0)= r [uty ~ > Py end =], 


ma 
‘ ge 
c,(0) = 1, Co(0) = M20, 
rae We 93/2 P. , Sma 3/225 , 


14 





where Q is the quadrupole moment of the deuteron, a and HL 
are the proton and neutron magnetic moments respectively, 


fee) 


and Py = f w?(r)dr is the D-state probability of the 
eieeron. 

Relativistic corrections to the deuteron scattering 
problem have been calculated by Gross [6], and Casper and 
Gross [18] and Friar [8]. These corrections account for 
relativistic modifications of the deuteron wave functions 
and the nucleon current incorporating not only Lorentz 
contraction in two dimensions, but also the fact that the 
equal-time wave function in the rest frame is not the same 
as the equal time wave function in the moving frame. 

schumacher and Bethe [7] and Elias et al. [15] have shown 
meee Che relativistic corrections of Casper and Gross [6,18] 


Z 


to the deuteron form factors, to order mee, 
M 


2 
ac LOr (L- 3) on the right hand side of EowacLlon il=7- 


enter as a 


Friar has shown that the relativistic correction enters 


il 


fmeewO Ways; first, a factor (l+t) ~ should multiply the 


right hand side afeeimemegekors hi-/; second, CE must be eval- 
uated at an effective a = q*/(1+n) . A comparison of the 
results wating both formulations may be found in Section V. 
In Appendix D, it is shown that by expanding Friar's 
correction in a binomial series for low q*, Schumacher 's 


4 


correction differs by less than 4 parts in 10° for the 


highest q? investigated. Consequently, the change in Ce 1s 
negligible at low ao. The relativistic corrections of 


Casper and Gross [18] in the form developed by Schumacher 


mealeare used in this work. 


JES 





When Equation II-7 is used to analyze the deuteron 


experiments, a scaling law for G is used [3]. 


ac 


The quadrupole term is dropped since it is less than 5x10" 


in comparison to the total at the highest value of q? 
reached in this experiment [38]. Equation II-7 may be 
written in the form 
Zee 


- a 2 a 2 20 eal 
E 054 (8) {Se +28 Gp ,t [bea (an) = alte 2 


ao 
dQ 


Z rug [1+2 (1tn) tan = ] and noting that 


Defining Ca = 3 


ao@ 
aN 


a D 
0, (8) Ga - 


where Gy i iemEeotal tOrm tLactor GL the deuteron, then 


2 
Gt = (G ~. (II-8) 


d geass 


) Cr (1 + Cg) Q - 


ra 


2 2 
ED ED 


Additionally, by defining ma as the deuteron magnetic and 


M 
16, 


relativistic correction factor, the following is obtained: 
ao. il ( 1 \ 
M. 14+C, ny 2 


meom Equation I1-8, the deuteron electric form factor, Gp ' 





d 
2S obtained in terms of Ka : 
1S 
aie 
2 ee ye nes (f) Gas yee 7 
Gas Ga Ky ae ie - (G,, ce ) cae (II-9) 
a ae Oy, (8) is p n 


6 





tie tatio Of the proton and deuteron electric form fac- 


tors, is expressed in terms of Gp and Gp as 
n Pp 
2 2 2 

Sp, (da? (G. (q3) 1p Some cay: ; 

ee ee ee ee iO) 
G Ca G ies E 

E E 

e & pe 


fenpbaning BGuations Ti-6, and II1-9, with Equation II-10, 
the ratio in Equation II-10 can be written in terms of 


measured quantities only. 


_! 
eZ 





1 
Mae) + G. (a2)lc aot 7% doP 
ie d E d E —— — -— 
n 7 d AL eel 
= |K Ky ——— (II-11) 
G. (q2) eoarey lee Ge (oy 
EP M M 


It is important to note that qq and oe which are the 
four momenta transferred to the deuteron and proton respec- 
tively have different values. The incident energy was 
adjusted to make the final electron energies from both the 
deuteron and proton scatterings equal. fThis insures that 
the transmission efficiency of the spectrometer and the 
detector efficiencies were constant for both gases, but the 
four momentum Geensifers were Slightly different. The 
deVries b' fit [19] for the proton electric form factor was 


used to adjust for this slight difference in momentum 


@eansiers. 
ee (q4) 
2 Y p 
p 2 
p ie G,, (q__) 


p P Jdevries 


The ratio in Equation II-1ll may be written 


vey 








dia 2 
G. (a4) + G, (a4) ie) ® (Se aS) 
n ee ee M iG p 
Se (q4) “gz \ doP of (oe) KE qs (q4) 
p AQ p deVries 
; 2 
Solving this for G. (q3) yields: 
n 
ay 2 
R Au Sy OEY 
Me 2 exp 1 p BP 1 
G., (qq) = Gp (q 4) R p> 5 ae Leo) (II-13) 
n p Mott K Ga. (ce) E 
M E d : 
p devries 
Rexp is the experimental ratio of the deuteron cross 
section to the proton cross section, and Serer is the ratio 
of ov (8) to of (8). Note that of (8) and Ke, were calculated 


at the incident energy of the electron on the hydrogen, which 


was slightly higher than the energy of the electrons on the 


deuterium. 
Writing 
d 2 
SE 2 . Ky Sp ‘o 
_. 2 |e 5 se) ee Se ee Se 
G R p 2 ae 
E Mott Ky, |G, (qq) 
P p deVries 


Equation II-13 can be rewritten in the more convenient form, 


: Ge 

2 ant 2 
G, (qj) =4n—-— = - 1.0} GL (qa). (II-14) 
E, 4d ce Eee od 


p 


fiyechis form Gp (q4) includes the effects of the Gross rela- 
n 


Pavastic corrections to the deuteron form factor. 
B, NEUTRON ROOT MEAN SQUARED CHARGE RADIUS 


In the Breit or "brickwall" frame, the form factors 


Gp and Gy can be separated. The Breit frame is the 
n n 


18 





reference frame in which the incident and final electron 
energies are equal, and the initial and final momenta of 
the deuteron are equal but oppositely directed [20]. 


In this frame Gp (q7) is the Fourier transform of the 
n 


charge distribution [21,22] and can be written as 


Z - iq-r 3° 
G, (q°) =f 9 GOs? = alse 
igs O 


The exponential term may be expanded in a power series 
about a = 0 and the coefficients of q” are the oe moments 


of the charge distribution. For symmetric charge distribu- 


tions, only the even moments are non-zero. Gp is given by 
n 
4 4 
Ge >qg 
2 ] 2 2 CH 
—= =— — + e if ea, 
Sg ta ) ean”? 6 <“cy? J Po a ree Cle sy) 


where <i is the mean square charge radius of the charge 


distribution. The slope of G, asa Pn Cite MeO £ q? for low 
n . 
momentum transfers is related to the mean square neutron 


charge radius by 


de 





Tit. SEAPERIMENTAL EQUIPMENT AND PROCEDURE 


A. NAVAL POSTGRADUATE SCHOOL ELECTRON LINEAR ACCELERATOR 

The. 100 MeV NPS LINAC is an electron accelerator similar 
in design to the Stanford University Mark III linear accel- 
waereor described "by" emedorow et al. [28]. The thirty foot 
accelerator is powered by three Klystrons, each of which 
receives up to twenty-two megawatts of peak power. The 
beam energy may be varied from 30 to 105 MeV with an average 
maximum current of about IA. The machine is pulsed at a 
rate of 60 pulses per second, and each pulse is about one 
microsecond in duration. 

The beam energy is selected by the beam deflection 
magnets and energy defining slits. The field of these mag- 
nets 1S measured by a nuclear magnetic resonance (NMR) probe. 
With the aid of quadrupole magnets, the beam is- focused to 
meecam diameter of about 2 mm at the target. A detailed 
description of the NPS LINAC is given by Barnett and 


Cunneen [24]. 


Deol ECTROMETER AND COUNTING SYSTEM 

The spectrometer, described in detail by Oberdier [25], 
1s a movable 120°, 16 inch radius, double-focusing magnetic 
spectrometer able to detect electrons scattered at angles 
between 45° and 150° in steps of 15°. Ten NE102 plastic 
scintillation counters span part of the focal plane of the 


spectrometer, where each detector serves as the source of a 


20 





counting channel. The width of the counters allows a momen- 
pam spread, Ap/p, of about 0.3% for each channel. Two large 
detectors, spanning the same field as the ten smaller 
detectors, served as dual coincidence backing counters. 
Thus an electron passing through both backing detectors and 
a front detector will register as one count on the Jorway 
Model 1836 Dual Channel Scaler located in the control room. 
A triple coincidence pulse will register on a scaler only 
feet OCCUrsS during the time that the counter gating circuit, 
gated with the accelerator gun pulse, is open. This pre- 
vents background counts from registering between electron- 
beam pulses. 

Scalers are capable of recording periodic pulses up to 
150 MHz, but due to the pulsed nature of the beam, the 
count rate was held below 60 counts/sec on all runs to avoid 
large count rate losses. The count rate correction was 
found to be linear up to 100 counts/sec and is applied in 
the codes BERHFORM and TRADFORM. It was decided to use only 
counters 4 through 7 in order to consistently follow 
poinc's Pe meeeecciire omiecesthe data points accumulated by 
Topping were used in the present work. Sufficient data were 
fakem SUCh that the ratios found from channels 4 through 7 
were well defined. 

The gas target cell used was a two inch diameter steel 
cylinder with 1.0 mil Havar windows positioned to allow the 
incident electrons to enter and the scattered electrons to 


ort. Both the Ha and D, gases were maintained at liguid 


PaAl 





nitrogen temperature and about 150 psig. The gas target 


system is described in detail by Topping [5]. 


C. DATA ACCUMULATION 

Before entering into the actual procedure for obtaining 
data several terms must be defined. A "point" is the set 
of data obtained at one spectrometer setting for either 
the proton or deuteron. It consists of the number of 
counts registered in each of the ten counters, the integra- 
Epon Capacitance and voltage, machine and spectrometer 
energies and the time required to take the data. A "peak" 
is the collection of data associated with one of the ten 
Saamnels, e.g., "peak 4" refers to the data collected by 
Memmcer number 4. A “run” is the collection of all the 
points of all the peaks for both the deuterium and hydrogen 
merone Value of gq? and scattering angle. 


Topping [5] reported results of 21 runs from ni = 0205 


2 


ww 


2 HicmaeSuLes ObenlS runs ana the 12 runs 


mo to 0.35 fm 
made by the authors will be combined to extend the data set 
mo 8.05 fm *. The data collected by the authors were mainly 
obtained under the same conditions and same methods as used 
Bee Topping [5]. 

All data for each point were recorded automatically (see 
Section IV-A). A run consisted of between 50 and 75 points 
first taken with hydrogen in the target and then repeated 
with deuterium. 

The spectrometer was first "saturated" at a field value 


corresponding to an electron energy of 100 MeV before data 


22 





were taken on either H. Or D. to insure that the same 
Meeeteresis loop*’was followed on both spectra. The spectro- 
meter was then set such that the elastic scattering spectrum 
was well above channel 4. The energy selected by the 
spectrometer was stepped down 0.05 MeV between each point. 
This process was continued until the energy selected by 

the spectrometer was about two full widths at half maximum 
below peak 7. Both the H., and D. spectra were obtained in 
the same manner. 

Two additional bits of data were recorded by hand, the 
Secondary Emission Monitor (SEM) efficiency and the pressure 
in the gas target. The SEM efficiency was measured when 
each of the peaks registered a maximum number of counts; 
therefore, a total of eight values were taken for each run 
and they varied less than 0.1%. 

The above data were always taken under conditions that 
were aS Similar as possible for the deuteron and the proton. 
The authors made runs from q? = 0.4 fm * to 0.5 fm *. The 
time required to complete a run in this experiment was 
Mreater than for the experiment of Topping because of the 
decrease of the cross section with ae and the fact that 
back angles had to be used exclusively. Consequently H. and 
D. spectra were accumulated on two successive days. This 
was a deviation from Topping's procedure, but this method 
gave data that reproduced within statistics. 

The beam charge integration was always set such that the 


hydrogen peaks contained from 10,000 to 20,000 counts in a 


Single channel. 
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As in Ref. [5], this work eliminates the need for 
ee ’ dof do P 
eounter efficiencies by calculating a ratio of so [oo 


for each channel. 
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iVe7 DATA ANALYSTS 


DATA REDUCTION 

The data collected on a Model 33 teletype was transcribed 
onto computer data cards compatible for use on the IBM 360-C. 
The data cards were organized such that there were two cards 
for every point obtained at a given spectrometer setting. 
Migs first card contained the counts for each of the ten coun- 
ters, and the second contained the number of background 
counts, the integration time, incident electron energy, the 
voltage integration, the energy setting for the spectrometer, 
the capacitor setting, and the integrated beam charge as 
measured by the SFM. The hydrogen and deuterium data were 
treated separately with FORTRAN programs BERHFORM and 
TRADFORM respectively which are listed in Appendix E. These 
programs are rivised versions of TOP#FORM [5]. They are used 
primarily to calculate the cross section measured by channels 
4 through 7, but they also make count rate corrections and 
charge monitor normalizations. Charge monitor normalization 
Was required because the automatic termination of the inte- 
gration for each of the spectrometer settings did not always 
stop the counting system at the same value of heam charge. 


The correction was 


V 
N(normalized) = N (actual) si SPER ICINS (TV-1) 
actual 
where een iSsyene actual voltage at the end of integration 
and Vee 1S an arbitrary standard voltage used to nor- 


malize the hydrogen and deuterium counts. 
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These programs were also used to plot the normalized 
count rate corrected spectra on the CALCOMP plotters at the 
NPS Computer Center. For Cio iast Plotting run, Cemtain 
additional data, which are necessary to calculate cross 
sections, (such as the radiation corrections and densities) 
were given arbitrary values. The count-rate corrected and 
normalized counts for both deuterium and hydrogen peaks 
were plotted. Using these plots, the position of each peak 
and the line shape of the electron spectrum was visually 
determined. A smooth curve was fit through the points and 
those points that were affected by energy instabilities were 
appropriately adjusted. Typical plots of both raw data and 
smoothed curves are reproduced in Figures IV-l, IV-2, IV-3 
and IV-4. 

The radiation correction program TOPRADCR [5] requires 
the following inputs: 1) the peak energies; 2) the pressure 
in the gas target for hydrogen and deuterium; and 3) AE which 
is the integration interval from the peak energy to the 
final spectrometer setting at which data was taken. These 
data rece obtained from the before ment ned plots. TOPRADCR 
also used the composition and thickness of the gas target 
windows and the virial coefficients of hydrogen and deuterium 
as stored constants. TOPRADCR calculated the Schwinger 
radiation, bremsstrahlung and Landau ionization corrections 
along with the density of the gas for each channel (peak). 
The necessity for the radiation corrections will be discussed 
later in Section IV-C. Average background was determined 


from the spectrum above the elastic peak. 
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Tne COtalsradaation and ionization corrections, the 
gas density, the SEM efficiency and the background counts 
for each channel were typed on computer cards and were used 
as input data to BERHFORM and TRADFORM along with the 
appropriate data decks for hydrogen and deuterium. These 
programs then recalculated the elastic cross section and 
the statistical error in that cross section for each channel 
considered. The calculation of the cross section is 
G@iscussed in Section IV-B, and the method for calculating 
Meatistical error is derived in Section IV-D. 

Through the use of BERHFORM and TRADFORM, four cross 
sections (one for each channel) were obtained for each of 
the hydrogen and deuterium peaks. These cross sections 
were then put in the form of a ratio (deuterium cross 
section divided by the hydrogen cross section),and a 
weighted mean ratio, ar was calculated. The electric form 
maietor of the neutron, Ge , for each run was determined 
using R.° The values of ane DPrOrenmeleeerieG form Lactors 
and the proton Pcutenen Mott cross sections and magnetic 
corrections were obtained from programs TOPROTON and 


BOPDEUTR [5]. 


B. CALCULATION OF THE ELASTIC CROSS SECTION 


The experimental cross section for each channel is 


CG Ao SIOe Sliema oie b 2 Sl _ 
exp Wi. N, aw o™ ster ; (IV-2) 


Si 





where ae is the number of scattered electrons counted 
Mea carerocular ehannel; 


= is the number of incident electrons, 


N is the number of target atoms per cm“, 


AQ is the solid angle defined by the spectrometer 
Siacs, 


Kees Kye and K. are radiation and ionization cor- 


mections (see Section IV-C) , 


N 1S given by 


ot Ns 
Ny = 7 P (iV=2) 





where fe) is the density of the target (g/em*), 
t is the target thickness (cm), 
No 1s Avogadro's number, and 


A 1s the atomic weight of the target gas. 


N. ais determined by 


in 
1 oe (IV-4) 
in e ee - 
sem 
where Q is the’ total charge of incident electrons, 


t? coulombs), 


e is the electron charge (1.60210 x 10. 
V is the voltage of the integrator in volts, 
‘© is the capacitance of the integrator in farads, 


E em +S the SEM efficiency (about 5.8%). 


In this experiment a ratio of cross sections is measured, 
so that several of the parameters of Equations IV-2 and IV-3 
which are not absolutely determined, i.e., solid angle and 
target thickness, cancel with a subsequent elimination of 


systematic errors. 
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N is the total number of scattered electrons, 
se 


Area under scattering peak 


Noo = Ginscounis Mev) (Iv-5) 
OF 


where E is the energy width of an individual detector. 
Pnieorder to obtain Oe. the number of counts in a par- 
ticular channel must be integrated over the energy range 


Of interest. 


ge dL 
_ 1 ws 2 
Nees a eae N ewe wlaa(ES E441) /0-0034877 Ba, lve) 
where N. is the number of counts in the channel for the see 
spectrometer setting, and E. is the energy of the ee spec= 


trometer setting. 
SomotiIng Bouations IV-3, IV-4 and IV-6 in Equation IV-2, 


the experimental cross section may be written as 





n-l 
l ii 
Pe oS! 2, Ba! o.00084 877E, (1.0+0.0039 (n-5)) Ae a 
Oo ec re a 
39 CV ot NT 
Bee a7 AQ 
A (iy 
where the expression 
a te oe (1.0 + 0.0039 (n-5)) (IV-8) 


expresses the electron energy seen in counter n in terms of 

the spectrometer dial setting Ea. Note that counter 5 sees 

the electron energy read from the spectrometer control dial. 
Count rate corrections were applied to N; Dyeueing a 


linear correction empirically determined from a thick target 


3 





experiment. The empirical correction is 


NG (recorded) 


N. (count rate corrected) = T.0 -0.00097 N,/t, ' 


where N,/t,; is the average count rate. Background was deter- 
mined by averaging the counts for each peak above the start 
of the elastic peak. Charge integration normalization was 
done aS given in Equation IV-l. 

Cc. is defined as the counts in the see channel, corrected 
for background (B), count rate and charge integration 
See iization. 

aa = 
sloje 


Ss SS Aen N,7t = CSS) 


The final expression for the cross section in terms of 


defined quantities is 


il = oare: (E.-E 4) KKK. 
2 1 sit lo mOmOOS 4s 77k. /E.) 
cle, ae © 0 (IV-10) 
exp ( ey, arte No 7 
ff Som (— 2 AN 
3 


where AQ is the solid angle and has a value of 1.831 x 10. 


steradians. 


C. RADIATION CORRECTION CALCULATIONS 

Both the proton and the deuteron cross section were 
derived with the assumption of one virtual photon exchange. 
This process is represented in Figure IV-5. For the deuteron 


the assumption of the impulse approximation is also made, as 


34 





re 
7 a 


; 


Feynman Diagram of 
Rosenbluth Calculation 


meuicemy = 
Pp a 
f Cf 
fe 
: "g 
PD 
j RF 
Y nf 
: BOR 2 


| | 
Feynman Diagram of 
Impulse Approximation 


Figure IV-6 








represented in Figure IV-6. However, in the actual scatter- 
ing process, more processes than those illustrated in the 
two figures occur. The experimenter must correct for 
several competing processes in order to be able to compare 
experiment and theory. These competing processes are 

1. The emission of real photons and the reabsorption 
of virtual photons in the field of the scattering nucleus. 
The correction to the cross section for this process was 
first derived by Schwinger [26] and is therefore called the 
Schwinger Correction. 

2. Bremsstrahlung emission in the target before and 
after the nuclear scattering of interest. 

See hscitation and 10nization of material in the target. 

iyeonrder £O Obtain the cross section for a nuclear 
Scattering event one would have to continue down to zero 
energy so as to count all the electrons which have some 
probability of losing all their energy in the three processes 
previously discussed. In practice the measurement of the 
cross section is stopped at some AE below the elastic peak, 
SO it is necessary to include a correction to account for 
the undetected electrons with energy less than E- =u. 
Where Ee is the energy of the elastic peak and AE is the 
energy from the elastic peak to the energy corresponding to 
the last spectrometer setting, usually two to four full 


mreeaths at half maximum. 
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C. RADIATION CORRECTION CALCULATIONS 
i “oCchvwndes Connection 

The Schwinger correction was generally about 15% in 
this experiment and is the largest of the three corrections. 

The Schwinger correction accounts for processes 
which occur in addition to the one virtual photon exchange. 
Schwinger's original work was later improved by Tsai [27] 
and by Meister and Yennie [28] with the inclusion of the 
effects of proton recoil. The correction used in this work 
was evaluated numerically using the expression taken from a 


review article by Maximon [29]. 


L 
6 : (hae ee 2 2 
J So 20 eee G7. \ 2a | 13 ca 
K. = @ = ono) 2 ¢ ts) (en “7 L 1 aa 4B fon =5 ) 


where 


E. and E, are the electron energies before and 
dbwer scattering, 


ais the fine structure constant, 
m is the electron mass, 
gq? is the four momentum transfer squared, and 


* gn (1-t) 


L. (x) = =| = dt is the Spence Function. 
0 


m  Bremsstrahlung Correction 
When an electron passes through material, it will 
radiate energy through collisions with atomic electrons and 


nuclei other than the target nucleus of interest. This 
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energy loss which is referred to as bremsstrahlung has been 
studied by Bethe and Heitler [30]. The energy loss per 


distance of material traversed is given by 


where XS is the radiation length, which is the mean distance 
over which the»releetron has its energy reduced by a factor 
et l/e. 

Following Butcher and Messel [33] the radiation 


length for a multi-element target was calculated using 


-( ; Jl: “ier 
x Vile Gien PO : 1 s 
oO 3 1 1 1 x 

1 O 


where x" 1S given by 





nA. 

SS SS er ee 

S 4a N r* AO. 0) OO. Less gt 3) 
Oo -e : TL 


where Nis Avogadro's number, 
Yr ais the classical electron radius, 


th 


Oo. is the density in g/cm? of the 1 element, and 


A. iS in grams/mole of the ith element. 


Hofstadter obtained the bremsstrahlung correction 
factor by integrating the probability function P(a) given 


by Bethe and Heitler. 


 (t/kn?) 
P(e) = ——S—————__ , 
ett eric) 
where P(a) is the probability that an electron traversing 


t radiation lengths will have an energy loss less than a 


Site! 





times its initial energy (aside from the ionization loss and 
the scattering from the nucleus of interest). The brems- 


strahlung correction factor 15s 
es r mee 
ape nP | Le () | 


where T is the thickness of the target (this includes 
the two steel 1 mil windows as well as the gas) 
in g/cm2, and 


E is the geometric mean electron energy in the 
target. 


In this experiment the bremsstrahlung correction is 

about 2%. 
fee Lonization Correction 

Electrons also lose energy by excitation and ioniza- 
tion in the gas and in the gas target windows. Landau [34] 
calculated the average energy that a mono-energetic beam 
would lose, and the spread of the energy distribution of the 
beam after passing through a thin target. From-this distri- 
bution it is possible to calculate the fraction of electrons 
which are not eae because of energy degradation through 
oma Zation. 

The fraction, P(AE), of the beam one does not detect, 


Memeause of 10nization losses is 


E le 
2 Ine 
W (AE) Zz N (E) an/f Mek sake, 


O O 


where E. is the energy corresponding to the energy of the 
last spectrometer setting, 


es is the incident electron energy, and 
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N(E) is the number of electrons detected at energy E. 


Landau [34] derived an approximation for W(AF), as 


have others [35], which is used in this work: 
W(AE) = lfw , 
where the value of Was provided by Borsch-Supan [36] is 
a QnA +A 
o = a(. - a4), 


where A = yrl = -0.4228..., 


Y is Euler's constant, 


AX = AE/E - 0.05 
iy 2, 
E =m 54 (oe ee T MeV, and 
Soi 
a 
T = target thickness in See 
n. is the number of atoms of the it® element per 


molecule of the target, and dis and a eae 
the atomic number and atomic weight of the 
ith element. 


~ 


iijeseoOnsazeation correction is thus given by 


_ | -1 
K, = [1 - p(AE)]™~ . 


For both the bremsstrahlung and ionization correc- 
tions, the target was treated as a homogeneous molecular 


compound of steel and either H, or Dos whereas in actuality 


2 


the target consisted of H, or D, gas (2 inches thick) at 


Z 
77°K and 150 psig and two 1 mil steel windows. The molecu- 
lar formulas used to approximate the target were either 


Fe Ho¢ or Fe Doe depending on the gas used. 
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Pe ERROR ANALYSIS 

im Systematic Errors 

In this experiment the proton and deuteron cross 
sections were taken with the same geometry, current, and 
final electron energies. At most a time delay of twelve 
hours between the end of the hydrogen and the start of the 
deuterium run enabled the proton and deuteron data to be 
taken within a 48 hour time span. Consequently, in the 
ratio of the deuteron to proton cross sections, many of the 
systematic errors cancelled; e.g. solid angle, scattering 
angle, SEM efficiency, target thickness, and absolute 
spectrometer efficiency. 

A systematic error which would not cancel in this ratio 
experiment is that which arises from impurities in the 
hydrogen and deuterium gases. The gases used in this experi- 
ment were respectively 99.5% deuterium and 99.995% hydrogen, 
Bemaetermined by their supplier [37]. However, the equip- 
ment used to establish these purities was able to detect 
ery 0.3% or more contamination of D. nd H. or H. in D.- 
The data taken in this experiment used the same bottle of 
deuterium and two different bottles of hydrogen. Consequently, 
the hydrogen data may have varying amounts of possible con- 
taminants depending on which bottle was used. However, if 
there were a difference in the contamination of He on Dos 
this difference would have caused a non zero intercept of 
mr versus q*. This difference was not evident 
as the value of G (q? 


EB 
n 


accuracy for all deuteron models tested. 


the plot of G 


— “== ZerOutco wrelin Statistical 
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2. Random Errors 

The statistical error is due to the Poisson distri- 
bution of the number of electrons scattered for a given 
number of incident See ee This error was calculated by 
a subroutine in the BERHFORM and TRADFORM programs for each 
individual peak and included the statistical errors made in 
determining the background subtraction. 

If X is the cross section, the following may be 
written 


M 
Phe—r (Constant) - 2 <{(N.=—- B) , 
me 


where Ny 1s the total counts (count rate and voltage 
corrected) in the kth Spectrometer setting, 


Bis the average background to be subtracted and 
1s assumed to be the same for each spec- 
trometer setting, and | 


Mis the number of spectrometer settings. 


fmeerollows that if N = Ny where N is the total number of 


counts, then: 
X - (Constant) - (N - M-B) 


mae Variance of the cross section is 


Z 2 
2 ox 2 9X 2 
ox = ON Nn” 38 [°° 
we = ieeneeany (N + M*o*) ' 


where of is the variance in the average background, and is 


equal to Bs Let Bin = Bs (where Bs 1s equal to the back- 


ground at the ith spectrometer setting) and let J be the 
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number of spectrometer settings above the peak at which 


J 
y Be 
| rei = 
background was taken. Since B = aa ae then 
J J 
On op Eon = ny By Sp Bp Fe 
J~ i=l an J~ i=l J 
maech in turn implies that 
of = (eons anu (N + nC =) P 


mm@meeene fLractional standard deviation of the cross section is 


1 
2 


x + Me al 
x 


x (N-MB) 7 

Besides the statistical (Poisson) counting errors, the 
random errors in this experiment originate from uncertainties 
in determining the density of gas in the target, variations 
in SEM efficiency during a run, energy shifts of the inci- 
dent beam, leakage currents in the beam current integrators, 
and small errors in the radiation corrections resulting from 
errors in determining AE. 

Uncertainties in determining the density of the gas in 
the target were caused by the variation of pressure in the 
target, that is, the pressure dropped while the liquid 
hitrogen was being pumped and then rose again and stayed 
Steady until the next nitrogen filling. These variations 
_ were almost identical in the hydrogen and deuterium runs 
and were about 0.3% of the measured pressure. The average 


pressure for each peak in the hydrogen and deuterium runs 
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was determined separately by either taking the pressure at 
~me start of the nitrogen ca} or at the end of the fill 
depending on which pressure was more constant. If the 
pressure at the start (completion) of the nitrogen fill was 
used for hydrogen, then the pressure at the start (comple- 
tion) of the nitrogen fill was used for deuterium. The 
errors in the measured cross sections due to the uncertainty 
in determining the density of the gas in the target varied 
from run to run, but were always less than 0.1%. 

The efficiency of the SEM was measured four times during 
one data run on each gas, one on each of the four elastic 
peaks. On each peak, five or more measurements were made. 
The fractional standard deviation of a single measurement 
was usually less than 0.2%, giving a standard deviation of 
the mean of the SEM readings of less than 0.09% for the four 
readings. 

Due to variations in line voltage, the beam energy was 
not constant, and as a result the scattered energy also 
changed slightly from time to time throughout a run. Also 
there was a slight instability in the spectrometer power 
Supply. The instabilities had the effect of shifting the 
counts from one channel to another on the counting ladder. 
This variation in the data was most noticeable on the steep 
Slopes of the elastic peaks where a low count in one channel 
waS accompanied by a high count in the adjacent channel. 

It was found that by closing the energy defining slits more, 


PiesVariation in the counts due to shifts in energy was 
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almost negligible. At the same time, it was noticed that 
the more the energy defining slits were closed the higher 
the background became. Thus, the ratio of real counts to 
background counts became less, which in turn increased the 
errors on the measured cross section. The position for the 
energy defining slits was adjusted to give the minimum 
combined contribution to the error in the cross section 
from the effects of energy shifts and increased background. 
Even after this was done, in some cases there still existed 
a need to adjust for obvious energy shifts. This was 
accomplished by smoothing the data points to the general 
elastic peak line shape. It is estimated that in the 

worst case these energy shifts could contribute an uncer- 
tainty in the average measured cross section of as much as 
0.2%, which the usual error in the measured cross section 
due to these energy shifts was ahout 0.1%. 

In the circuits used to measure and integrate the beam 
current there existed slight leakage currents. If these 
leakage currents — NMOowmecoOnseamt threugqnout a run, then 
they would introduce errors. These currents were measured 
and found to be constant as well as negligible [5]. Asa 
result the only error attributable to these leakage currents 
existed when the electron beam current varied. The electron 
current was usually held constant to within 20%. This 20% 
variation would cause an error in the ratio of the measured 


cross sections of about 0.023%. 
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Pothier Seurce Of uncertainty 1s the inability to 
determine the exact position of the peak of the elastic 
spectrum from the plotted data. This uncertainty is 
usually about 0.01 MeV, and causes an error in the Schwinger, 
bremsstrahlung, and ionization corrections. The error in 
the measured cross section due to this AE uncertainty for 
all three radiation corrections was of the order of 0.05%. 
An uncertainty in the final energy at which one terminates 
the integration of the elastic spectrum, contributes approx- 
imately a 0.03% error to the measurement of the cross 
section. 

A list of all the errors discussed can be found in 
Meere IV-l. Excluding statistical error (since it varies 
Soemuch from run to run), the average error in the cross 


section ratio, S is about 0.24%. The total 


Avg. Others’ 


error of the weighted mean ratio R was then formed from 


2 mee: 2 
eal 7 Duet a3 Pav. Others (Iv-11) 


One significant source of error has been reduced since 
Hepping [5] performed his ratio experiment. It was found 
that if quadrupoles 3 and 4 were used instead of 1, 2, 3 and 
4 that the background was effectively 50% lower. The most 
plausible explanation for this reduction of background is 
that quadrupoles 1 and 2 defocused the electron beam, caus- 
ing the beam to hit the beam pipe and thereby increased the 
background. The data in this experiment was collected with 


only quadrupoles 3 and 4 on, and as a result, the peak to 
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vac, Vo 


RANDOM AND SYSTEMATIC ERRORS IN CROSS SECTION DATA 


—————— 


Errors Average Values of Errors 


Bearcistical 0.3% 
Gas Density 02 L% 

perm Efficiency 07093 

Energy Shifts 0.1% 

Leakage Currents 0.014% 

E Measurement C205 s 


Peeaiation Corrections) 


E Measurement 0.033% 
(Integration) 





Average Error 0.24.3 
(excluding counting statistics) 


background ratio of counts was much higher than in the 
previous work [5]. 


The uncertainty of Gp arises from the errors made in 
n 
measuring the individual cross sections. Let R= 54/5), 


where o, and 0, are the experimental cross sections for the 


d 


proton and deuteron respectively. Let oe and S. be the 


d 


standard deviation of the proton and deuteron cross sections 


respectively. Gp is determined from 
n 


= ae 2 = 
os (6, icy, ) a - 1.0} G. (qa) . GEV = 152)) 
n d p E p 


and Gp /Gy can be determined from 
d “Pp 
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Pp Geecs (c=) 
Gc. /G Sah Omott Sage Hf Pp 
E E d KP fe 


5 (IV-13) 
d le Omott M E,, ‘Ip’ 


Solving for the standard deviation in R, denoted by Spe 


memento the errors in both cross sections, one obtains: 


2 ah iy) Wl _P s 
Sp = R 5 + -P , (IV-14) 


The standard deviation of the ratio G., -G, , denoted by 
cl 2 

Se /G, is obtained from Equation IV-14, 

d 


p 4 Cree 

0 Ine E p S 

Sc /G,, = eee = ——— R , (IV-15) 
d Ta 1 Sp ‘Ip’ 2VR 





From Equation IV-12, one can determine the standard 


deviation in G. , Genoted by Sc to be: 
n E 
n 


G (qa) G 
S eS : (IV-16) 


Os 
v8) 


mumorder to combine the results obtained from four different 


channels in one run, let R. Gen@eorvenec Vaclo Of Cross sec-— 


th channel and let S denote the standard 


mens in the i R 
ch 


1 
Bemmation of the ratio in the i channel. It follows that 


the weighted mean ratio Re can be calculated from 


7 2 

_ Ri/°R, 

R = =———_——_ ,, (IV-17) 

m fi 
; [1/s 

jean 

=4 1 











48 





where the standard deviation of this weighted mean ratio is 


calculated from 





I 
SR - ee . (lve Iss) 
» [1/s 
: ie 
i1=4 1 
R, and Sp are the experimental values of the ratio and 
m 


@: the standard «deviation to that ratio. 
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V. EXPERIMENTAL RESULTS 


A. NEUTRON CHARGE FORM FACTORS 

mie 2) runs by Topping [5] and the twelve runs by the 
authors are reported here. Each run consisted of four 
measurements of R, the ratio of the elastic electron- 
deuteron cross section to the elastic electron-proton cross 
section. Of these 132 measurements of R, 9 were not 
included due to reasons given in Ref. [5]. 

Table V-1l gives the q? transferred to the deuteron, 
the scattering angle, the experimental ratio Ru! and the 
standard deviation of as determined both for the total 
meme and the statistical (Poisson) error for a given run. 
Also included are quantities that were determined by cal- 
culations discussed in Section II. The ratio of the 
Pem@eeton electric form factor to the proton electric form 
mieeor, derived using relativistic corrections by both 
Friar and Gross, Sys also listed in Table V-l. 

Tables V-2,3,4 and 5 present the neutron charge form 
factors SHeeRCE using the deuteron models LF1,5, and 15 
with the relativistic corrections of Gross and Friar. The 
deuteron structure factors used are presented in Appendix D. 
The best estimates of the error in each of these 33 values 


of G,. for each LF model and relativistic correction is the 
n 


error obtained using §S from Equation IV-ll. An average 


total 


; 2 
G 
E and the associated Son at each q° were calculated 


using Equations V-1l and V-2 respectively. 
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A plot of the average Gy, (with its S ) versus q? isp Ee 


n 
muaed for each of three LF models and both relativistic 


one 


@orrections. 


All values of Gy, were calculated on a Wang Model 700-C 
n 
meeeulator using programs written by the authors. The best 
versus Ge a/aq? G 
lg : E 
n n 


obtained from a linear fit weighted by Sot for all values 


estimate of the slope of G , was 


St G. - For comparison purposes the linear fit was also 
n 


memeed through zero for each plot. 


The average value of Gr aie Cac q*, Useag Lomon- 
n 
Feshbach Model #1 is given in Table V-2. 


The average value of G at each aa using Lomon- 


EN 


Feshbach Model #5 is given in Table V-3. 


The average value of G at each a. using Lomon- 


E 
n 


Feshbach Model #15 is given in Table V-4. 


The value of Ga calculated from Lomon-Feshbach Model 
n 
#15 and its =e for each of the 33 runs is given in 
n 


Table V-5. 


of 





i. ROSENBLUTH PLOT 

At q? = 0.4 fm 7 deuteron and hydrogen data were taken 
at 90° and 120° so that a Rosenbluth plot could be con- 
Structed. The theoretical interpretation of the Rosenbluth 
plot is discussed in Appendix A. Table V-6 gives the 
experimental data and derived quantities used in the Rosen- 


bluth plot given in Figure V-7. The use of the scaling 


law 1s compatible with the data. 


~~ 


— 


See CHARGE RADIUS OF THE DEUTERON 

Using Method II given by Schumacher and Bethe [7], the 
deuteron structure rms radius was calculated. A more 
detailed explanation is given in Appendix B. The structure 
radius of the deuteron for LFl, LF5, and LF15 was provided 
by Schumacher and Bethe [7] and was used in a plot of the 
2 versus 2 (deuteron structure 


_ 
dq? EnJq2=0 Ea 
radius). This plot is Figure V-8 and as predicted by 


experimental 


Schumacher such a plot gives a straight line for all models 


with the correct binding energy. The intersection of this 


mime with the value of |, G 


dq” ae 

accurate neutron-electron scattering work yields ce . The 
d 

value for the deuteron structure rms radius was determined 


determined by the very 


from Figure V-8 to be rn = I 966s. + 020045. 


a2 





is) ellu, Gy Wacoal & 


EXPER IMENIAD RATIOS AND GRESULTS OF 
MODEL INDEPENDENT CALCULATIONS 
















Exp. hat1o 


4 Sra Deve 5 
_ ex om R Ge (ae 
Ree Iieyegsill 2kgte p 
P (stat.) 


exp 


~9654 0.3026910.99400;0.99404 


(WoOOO) IO. SLs 
9724 0.0038 |0.16837/0.99396/0.99404 
(Ono OS Ol Onli G7c 
-9390 0.0042 |0.81525!0.98789}0.98815 
(Or00S 5) a) Orciog 57 
.9331 0.0034 |{0.81525|0.98789}0.98815 
(020024) GOe sie Sy 


0.0046 0.81525! 0.98789!/0.98815 


9262 
86757] 





















2 d (Friar) 
qq ieee ss) G Total 
E 
a 2 ae, oe 
(fm ~) Ka (Friar) Gp —(stat) 


Error 


OFMeieco 0. UUme) 0.965510. 002s 
(0.0020) (0.0020) 


0.96934+0.0019!0.9693+0.0019 
(0.0015) (Oneo Oulkes. 


0294752050071) 0.947520.0021 
COP OOnrs ) (Ol OOiS) 


0944 02 0087 10.944520.0017 
(Oe 00A2) (G00) 2) 


Cane O002310.941020.0023 
1c LE) CORON, 


5) 6) 





TABLE V-l (Continued) 











Ee. Ratro 






















Std. Dev. 


d mott 
O of R 2 2 
R= —P |rotal and e Gp (4) | Gp (4@) 
oP (Stat. } Alc P P 
MOoce 


0.0049 Or 6152510... 98189 WW. 28S15 
CO 0031). 8677 


0.0051 10.70815 1N.97595| 0.97653 
On .MiNIaSYelN .FE120 


() 01017 6 O. 7OeH oa .97595 Oh, 27 Gira 
CORO 7 2 BO. el 9 


0.0049 SOAs) . 7 58 2 0.97633 
(OmO042 ) Wl t2650 


0.0049 Nee 2434) ..97569 Gees 3) 
(9.90042) (0.24871 


O.20i052 re 2a) 2.9) 569 NE SPAN Gy 8) Ss 
COOa6 ) Th. 249 7a 







a4 (Gross) Total 
E 

_9 d and 

(em ~) (Friar) r ~ (Stat) 








0. 93892020070 
(O02 OGG ) 


Ono 68 920002 0 
(0570 0nE Gs) 













Pood o=(). OZ 8 
(0.0024) 





5884620 70023 
(0. 0024) 












0.888720.9041 
0 0032) 


~Oeee 20s M4 | 
Coins ah) 








peOZ k=) 026 
(ORO O22) 


Arse) 210 selG (QO) 2S 
(0.0022) 















Orc 10 20026 
MOR O0IZ 2) 


8S (Oro 0 26 
CO O22 ) 








Memeo St. 0028 
(6.0024) 


sooo S20. 0028 
Corie) 










54 





iabbe Vol (Coatinmed) 









EXD. Ratio 






Std Dev. 
On TK 
Total and 
(Stat. ) 













SO) Lars. 
2CEOG) 








OO) 2 
. 010199) 








O82 
ay 9) 








.0056 
HO 52 ) 











0.0049 
(0.0043) 











0. O09 
(020052) 











(Gross) 
G Total 
—d , and 
G ~ (Stat) 
EEror 



















Seen OGG 5 
(0.00 s18) 


s COZ OO O01 = 0). 0 CG 
(Ozi0 Osis) 









POO Ss 0055 
(0.0054) 


~2901820.0055 
(0.0054) 











soe vo O04 S 
(0.0944) 


.88962£0.0N45 
(0.0044) 








peg = (igs 
(07,0028) 


cos O00 2i0 
Or, OZ ey 














.842620 0030 
COMGVay 


.8426=0.0030 
(0 ..NN26) 
















.853820.0031 
COL 020) 


852 C= Or OO: 
Cea») 





ID 





TABLE V-1 (Continued) 


ae 
_—_ 


0. 1433670 .96S6 5 








Eee. Rat1o 





Std. Dev. 
@f R 
E@itad.) ang 
(Stat. ) 













Ns Jos 6 







BW oy a LE 


<1 59e5 







[93055 









.95400 






.95400 







2 (Gross) 
da 
an *) (Friar) 











Hec42320 70064 
kOx 


Oe tss 0.00 s4 
(ORCn3:L) 









.99845 
99844 











~8449240, : 
MORO 3 7) 









1.844920 .0034 
CORR 7) 


Be) Soe. 
.99844 













Gm3iz4220 0081 
(0, 0027) 






08624250 .0032 
KO O02 7) 


21S Ste Ae) 
[29818 


















Wied ola) 
LOR C032) 


Nos 261 08 5 
(0. 0022) 


pee Slee alee 
SN) sides 
















feos 7+ 00033 
oe 


().. ol 7 2001033 
(Oe Os 0) 


2989 33 
SNS) IS) dts 






















a0 6 + 0. 
(0. 


0. c0R6 2020034 
(Cro ter!) 


omnes: 3 
02997918 





5,16) 





2 

dq 

(fm *) 

0.40 |120° 
120° 
120° 

0.477 | 90° 
90° 
90° 


OF 


0 


OF 
OF 


OF 


TABBE V-1 


Exp. Ratio 


d 


R = —2SP 


(Gross) 


(Friar) 


SSNS 
98936 


Seekers: 
7920)36 


5 oles 8) burs! 
293) 36 


O97 4 


So) oe 


99754 
SEES 2 


99754 


Std. Dev. 

of R 
Tonal and 
(Stat. ) 















219004 


Ses) 


oo a 


2389 35 


Si So 


Mrehove).2) 3 


(Comt imaued ) 













OF 


a7 


Geez ce 9 
Oe Beavess) 


ONS Fare 
01739462 


a3 2c0 4 
Oe Ser ge 


On coe 2 
il O06 2 


Clete chs y5 20”. 
1.09679 


OS 552 2 
1 GOS Fe, 


(Gross) 
Total 
d F and 
~ (Stat) 


BrLror 


Bis 020022 
(0.0036) 


780 200039 


0G. 010316) 


esos 9 


(0.00 3.68 


~/ (3900s | 


PUM ee |) eal 


(OS 002% 


718220 .. 0By, 


CORO 0 S29) 


5 to) IAB Syanl@) 
(0. 


wo Oa 
(0. 


mei) SO) 
(0. 


Ee ale 
coe 


. ho Ome 
(Cx 


mee 2 ee, 
COs 












































TABLE V-1l (Continued) 











Bre. Ratio Sedna 


of R 
Total “and 
(Stat...) 












(Gross) 
Total 
d a and 
(Stat) 
Brror 








TST ee 
(n. 


seu 0) aCe 6 
Con W133) 





















Pawo tO, 
(Of 


7635 0 O0S6 
Mae, O10 33) 









.768540.0037 
(0.0034) 


pe Gave = 0) 4 
(al. 













wo6OS0: 
(Or 


RT lStens aly, (lg exs 
(0 .. 08 59) 


2y8) 





TABLE V-2 


AVERAGE Gp AT EACH q? USING LOMON-FESHBACH MODEL #1 


Nn 


Gross me relativistic Friar's Relativistic 


q? Correction Gerrection 
(fm “) Cy Soe Cy Som 
we. 05 .0002 nOOds .0002 a ONC lies) 
oe 10 nO? OO ual OWE 7 SOOnE 
0.20 O27 001s - O27 5 ova 
or 25 70070 .0034 .0069 -0034 
0.30 SOO76 .0019 C075 CIC isle!) 
0.35 SoCo OO 7 .0092 O02 
0.40 -0084 . 0020 « OOS .0020 
Oe /7 7O098 /0025 20095 - 0025 
a 50 FOLO6 OO {OIL .0023 


A linear fit to the data in Table V-2 gives: 


G. = (0.0003 + 0.0010) + (0.0208 + 0.0040)q” 
Ber Gross’ correction and 
Gh = (0.0004 + 0.0010) + (0.0200 : 0.0040) q* 


foeenielai S COrrection. 


A linear fit of the same data forced through zero gives: 


GOGO 21 oes 0.0022)q" for Gross' correction and 


O 
tl 
io 
T 


COO 2 13 0.0022)q° Police Secor rection. 
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Rapin. Vs 


Hiren's Relativistic 


AVERAGE Gp AT EACH q* USING LOMON-FESHBACH MODEL #5 
n 
9 Gross' Relativistic 
gq Correction 

(fm *) S Sg, 
n n 
O05 0.0002 O20 
0.10 Oro 7 Ooo Tal 
oF 20 0.0028 Oreo 1k 
O25 OF 01017 0 0.0034 
wn 0 0.0076 0.0019 
i 35 0.0094 Oe ON BY 
0.40 0.0084 0.0020 
0.477 O00 99 0.0025 
0 OOO 7 0.0023 


A linear fit to the data in Table V-3 g 


G. (a7) = (0.0003 + 0.0010) + (0.02 
n 
for Gross' correction and 
G, (q*) = (0.0004 + 0.0010) + (0.02 
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heme wal Seeerrect om. 
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A linear fit of the same data forced through zero gives: 
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bl 
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TABLE V-4 


AVERAGE Gy AEE ACH q* USING LOMON-FESHBACH MODEL #15 
n 
5 Gross' Relativistic Pebaees Relativistic 
q COGBeSerLoOn Correction 
—-? ee aaa cr 0 SS a. 
En En Ea En 
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A linear fit to the data in Table V-4 gives: 


Gy (q*) —meoeOUO2 = 0. 0010)) + (OMO01l77 = 0.0040)q* 
ms for Gross’ correction amd 

CG. (q7) = (0.0003 + 0.0010) + (0.0170 + 0.0040) q° 
E 


EGG ELlar S COLrrecti1on. 


A linear fit of the same data forced through zero gives: 


a 
iQ 
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VOniOMES 5.52 0.0022)q for Gross' correction and 


(OZ 0180). + 0.0022)iq" EOL etiam es correction. 
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TeELEF Va5 


ie AND at 


CALCULATED FROM LOMON-FESHBACH MODEL #15 


FOR EACH OF THE 30 RUNS 


5 Gross' Relativistic Hewat: GeRelativistic 
gq Correction CGEREECtion 
-2 
En 5), Bn En 
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Var CONCLUSIONS 


mea ream tC OltLmehembles measurements of the ratio of the 
elastic electron-deuteron cross section to the elastic 


electron-proton cross section, nine values of Gp at dif- 
n 
ferent values of q? were determined for each deuteron wave 
function considered. From these values of Gr and aa a 
n 


value of d was determined. 


| 
eens) 2= 9 


Relativistic corrections proposed by both Friar [8] and 








Gross [6] were investigated. It was found that, for a given 

qeuteron wave function, = 5 Gy 5 was smaller using the 
dq n'iq =0 

fame COrrection. Although this pattexn was consistent, it 


was also found that the difference in the two slopes was 
always well within one standard deviation. No attempt was 
made to select one form of the relativistic correction as 
best because the difference was so small. 

All wave functions considered (LFl1, 5, 15) were in 
agreement with the very accurate thermal neutron slope of 


7 The values of es G 


dq 


e093 + .0004) fm E 2_, obtained 
njiq =0 


meamg Gross’ relativistic corrections were: for LF1 (0.0218 + 


Mme 22), £Oor LFS (0.0220 + 0.0022) and for LF15 (0.0185 


I+ 


0.0022). The values obtained using Friar's corrections 
fete fOr LF (020213 + 0.0022), for LF5 (0.0214 + 0.0022), 
emo tor LFI5 (0.01380 + 0.0022). LF15 was found to be in 


closest agreement with the thermal neutron slope. 


2 





Pa tom Gremio ssireuecture radius of the deuteron was 
also obtained from this data. This value was determined 


by comparing the slopes given by each wave function. The 


structure radius was determined to be Ce (1.9665 + 0.0045) £m. 
A Rosenbluth plot was made at q? = 0.4 ee The results 
of this plot are discussed in Appendix A and show agreement 
with the use of the conventional deuteron scaling law. 
If this experiment were to be repeated or extended it is 


recommended that data at higher q? be taken in order to 


determine more accurately the coefficient of the q’ term in 


the expansion of Gp meno enarx Cc a COeLiicient for the 
n 
qa: term was taken from work done at DESY [38] and it was 
femnd that at q? = 0.5 fm“ this term decreases Gy, by 
“0 


about 10%. The effect on the slope of the Gp curve is also 
n 


Smown to be about a 7% increase. 
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APPENDIX A 


io eROmeN Biwi POT 


Based on the one photon exchange, the cross section for 
electron-proton scattering was derived by Rosenbluth [10]. 
Using the same assumption (one photon exchange) Jankus [14] 
calculated the deuteron cross section. Jankus' result is 


given in Equation II-7 and is reproduced below 


aa aia 2 2 


8 Z 28 
as? hy, oad Eq g O 3 Ms 2 


where ee can be written as 


G2 = em (q*) + 8 nc G2 - 2 n G2 1+ iaeeene 
Eg 9 3 Ma 


No] D 
Sas 
Soe ee 


Semce the deuteron cross section can be measured at different 


scattering angles while at the same four momentum transfer, 


one can see that a plot Oye er versus (1 + D (7 hea teaiee >) 
should give a straight line with 2 n a as the slope and 
: d 


3 

Z BO 2 eer 
Gp + 9 (7 Go 2omelemei intercept me hy /PplOttang the deuteron 
d d | 


meocal form factor a versus eer Ci anban a], the adequacy 


of the conventional scaling law can be tested. Such a plot 


has been made for q? = 0.4 fm? 


anoeewean bee Lound ean 
Section V. 
SMe cCann@emise soquacion LI-/ directly since it contains 


the absolute cross section of the deuteron. However the 


deuteron cross section can be obtained by using the hydrogen 
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cross section along with the deVries b' fit. The deVries b' 
cit [19] is used because it is the best fit to the absolute 
proton cross section. It is used as a normalization factor 
with the experimental cross section obtained from the data. 
timat is, 

ao? doP 


= deVries exp 


Therefore, 


doP doP 


oH deVries an exp 


N is a normalization factor for a given experiment at a 
given angle and energy. Since both the proton cross section 
and the deuteron cross section were measured under the same 


experimental conditions, the deuteron cross section is 


do® - dot 
AR 7 AQ 
exp 
This implies that 
ao" | _ | /aoP ao? | [aot 
AQ AQ AL? \d 
deVries EXP |\ ePs)e) 


Now putting this in terms of the experimental ratio R, 


where R is given by 


R = dot doP 
AQ AN g 
exp exp 


@me Obtains 


cS 











dct =a doP 
aN oe deVries 
and recalling that 
dot = dot GZ 
dQ dQ ote d 
one gets the final expression for ce as 
ee = doP do R 
< oY deVries ae mott . 
doP 
where aq” was calculated from the BERROTON program 
deVries 
dot 
and aq” was calculated from the BERDEUTR program and 
mott 





Emme ratio Be is the weighted mean ratio given by Equation 





ay-17). 
ites statistical error in G4 was calculated using the 
following equation; 
jac 
S — a deVries = 
Ga do m 
dQ jmott 








where Sp 1S given by Equation IV-18. 
m 
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APPENDIX B 


RMS STRUCTURE RADIUS OF DEUTERON 


Pyestiostituting the deuteron charge distribution for the 
neutron charge distribution, it can be shown that an expres- 


sion similar to Equation II-15 for the deuteron is 


applicable. 
2 4 
<Ye> <> 
ae qd 2 a4 7 


where ae is the mean square radius of the deuteron charge 


Meter OuUtion. in principle, the deuteron form factor could 


Memrit to this polynomial form and the slope of G,(a°) 
ire oe 





versus q* in the limit as q? approaches zero would be - 
However, for the deuteron, the coefficient of the ae Sian 


not well determined and the choice of <r4> eacatlyeatiects 


d 
the value of <ro>. In order to determine re data at 
higher momentum transfer is necessary and its determination 
| 6 


requires knowledge of the coefficient of the q term. 
Fortunately, Schumacher and Bethe [7] have suggested a 
procedure to avoid this problem. The method consists of 
using the experimental data and the function cco) associ- 
ated with a deuteron wave function to calculate A in the 


expansion 


6. (66) = (B-2) 
ED q q 


Wii 


oe 





i the deuteron structure factor Ch is also expanded for 


small sa the expansion is: 


2 2 1 4 4 
Srp? q- + 130 <ry? 8 (es cee (B-3) 





fr 


2 _ = 
C.(q°) = 1 


where eae may be defined as the mean square radius of the 
Geuteron structure factor. 
Schumacher and Bethe observed that if the same data set 


memused to Calculate A from various C.. functions, which have 


2 
E 


E 
meeererent values of <ri>, the plot of A as a function of 
<ro> forms a straight line, if one uses only Ch functions 
corresponding to deuteron models which have the correct 
binding energy. Then, if one requires that A have the value 


deduced from the neutron-electron scattering experiment, the 


value of <r2> may be determined, as shown in Figure V-8. 


E 
The error in <ré> is fixed by the intersection of the pro- 
jection of the errors of A, and the slope obtained using 
2 


models considered in this experiment. The value of <r 
is given in Section V. 
Finally, one may note that the mean square charge radius, 
2 


<r4> may be found from: 


ee 2 2, 
5 = <r? + ee Oy ae 
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APPENDIX C 
A 


gq DEPENDENCE OF Gp 
n 
In Chapter IV, the slope of Gp was determined assuming 
n 
eeiinear fit of Gy, versus q’, but there was some evidence 
n 


Mmiat the inclusion of a quadratic term may be more appro- 
priate. Therefore, a fit with both linear and quadratic 
terms was made to the data, but it was discovered that 
mem a fit produced a quadratic coefficient of little 
confidence. 

Since the value of the coefficient for the linear term 
is not strongly dependent on the value of the quadratic 
Beetticient, for this data, a value for the quadratic 
coefficient was assumed using the work of Galster et al. 
fees Galster et al. obtained the value B = -0.0036 fm" 
for the quadratic coefficient and this value was used in 


application to this work's data. A recalculation for the 


quadratic coefficient of the analytic expression, 


Gp = Ag? 1S Bq? - 
n 
emews that the quadratic term contributes 5% to G, at 
n 
a eo. 355 fm Z and 10% at q? = 0.5 fm a With the assumed 


m@uadratic coefficient, the slope of Ga increased hy 73% 


n 
from that obtained without the inclusion of the quadratic 
Perm. Also, the deuteron structure radius is decreased by 


0.15%. Note that the changes in the slope and the deuteron 
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Piructtire radius are within the statistical accuracy of 
the value previously reported in Chapter IV. A comparative 
feote Of values utilized in the linear and quadratic fits 


can be found in Tables C-l and C-2. 
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COMPARISON OF G 


ae 
(fm a “EL 
ie OS 0.000 
Vinee Oru 
ie 20 Oen0 012 
i. 25 0.007 
o. 30 OnOg7 
o.35 Oe (Cs 
0.40 0.008 
OF 4/7 0.009 
te. 0 0.010 


For LF1l assuming B 


A 
Hem LFS assuming B 

A 
For LF15 assuming B 

A 


Pesos Cal 


4 
RE AND Gp = © | 
n n 
LE 5 LFS 
4 4 
Gy ao | Gp Gp -Bg 
n n n 
0.0002 G2 oc ged eo 0i18)6) 
00037 One0ss 02.0033 
0.0029 G00 ZU 0.0021 
0.0072 C2006 L 0.0063 
0.0079 0.0066 0.0069 
0.0098 0.0082 0.0086 
0.0089 Owen, O 0.0076 
O20107 0.0083 0.0091 
0.0116 0.0090 0.0098 
4 

=e) 36 52m 

GmoOz32 + 820022 

- 0.0036 fm* 

Gage ec et no 0022 

- 0.0036 £m4 

O2O19 8. = O2002 3 ‘ 
TABLE C-2 


COMPARISON OF THE DEUTERON STRUCTURE RADIUS OBTAINED 
PeOlvmian LINBAR AND QUADRATIC FITS TO 


G 


ASSUMING GROSS RELATIVISTIC CORRECTION 


E a 
n 
y) 2 4 
G = Ag G = Aq + bg 
Ey En 
<ro> BEe67 = 0-018 3.858 + 0.018 
1.9665 + 0.0045 1.9642 + 0.0045 
rms 
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APPENDIX D 


COMPARISON OF GROSS‘ AND FRIAR'S RELATIVISTIC CORRECTION 


Berti tne Gross and Friar relativistic corrections 
Wereroly the right hand side of Equation II-7. The Gross 


metativistic correction factor is given by: 


2 \2 2 4 
l= =, Se ee (D-1) 
8M 4M 64 M 
p p p 


emeePriar'’s correction factor is 


9 =. 
i+ —41, 
4M 
p 


Expanding Friar's correction in a binomial series to second 


2 
order in So , the following is obtained: 
4M 
p 
p a 2 4 - 
Lot thors si = et See, (D-2) 
4M 4M 16M 
pe p Pp 


The difference between CiewEWOmbeldelvVistile correction 
factors is approximately given by the difference between 
Peet iOn D—-2 and D-l. If this difference is given the 
symbol A then, 
q* i aq 3 q’ 
Az — - =—|r=- = ° ee) 
M 4 |16 64 64 4 
M 
p p 
The difference is seen to increase with increasing a 


Pet. q? = 0.5 fm“, which is the highest q? imvestigated, A 
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5) 


feos 10 with Friar's correction always being greater. 


2 2 


Pimce the corrections are both about 0.997 for q° =0.5 fm “, 
bhenethe diiference in these two corrections is about two 
fares in Or: Table D-1l gives a comparison of both correc- 
mons Lor all q*'s investigated 
Friar [8] also shows that Car the electric structure 
factor of the deuteron, should not be evaluated at q? but 
| 
at an effective gq given by 
5! 2 
es  . (D-4) 
1 + 4, 
4M 5 


This also creates a very small change as seen in Tables D-2, 
Deoeand D-4. Friar's first relativistic correction opposes 
the tendency of the use of the effective momentum transfer, 
so the net difference between the use of the Friar and Gross 
Berrections is smaller than four parts in low. Also 


included in these tables is a comparison of the slope of 


Gy versus q? using each correction. Again the difference 
n 


is seen to be negligible for this experiment. 
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Tass -D= i 


COMPARISON OF RELATIVISTIC CORRECTION FACTORS 


gq? PY jae --s Gross' ey 
Relativistic Relativistic 1- = 
-2? : Gross 
fm Correction Correction 
205 .999972 .999972 -.000000 
0.10 .999447 .999448 -.000001 
0.20 .998894 .998896 -.000002 
O25 .998618 .998621 -.000003 
0.30 .998341 -.998345 -.000004 
te 35 .998065 -.998070 -.000006 
0.40 -997788 .997796 -.000007 
0.477 .997360 .997370 -.000010 
0.50 oo 2 35 ~997247 -.000011 
TABLE D-2 
COMPARISON OF Cmca AND Ele or OR GEE 
' 
2 
Ca Clmiae) 
E 
aC C (q*) q? C (q? ) l SS 
E E 2 
Ce lq) 
05,05 .96882 0.0500 .96882 -.00000 
0.10 .93967 0.1000 93969 -.00001 
0.20 .88668 0.1999 .88673 -.00006 
0225 ood 7 0.2498 4 13\\6) P26 -.00009 
0.30 .83961 0.2998 . 83972 -.00013 
0.35 .81796 0.3497 - 81810 =r) (0)(3)0) ea 
0.40 .79742 0.3996 . 79760 -.00023 
0.477 . 76766 0.4765 . 76790 -.00030 
0.50 . 73929 0.4997 775955 -.00034 
Slope of GE, versus q? using Gross' correction is 
Oe0218 + 0.0022. 
Slope of Gp versus ae wlisuneemrrac's correction is 
Hoes + 0.0022. 
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TOS) 
70 
x20 
25 
30 
“0 ey 
0.40 
0.477 
i> 0 


oO oOo @ &€& yo 


y 
CoG ) 


96881 
se) 2216) [6) 
moe O65 
86244 
go59 56 
ore SU) 
ony 5 
ORO 
Sess) 7A Ml 


PABIGE. D— 3 


COMPARISON OF eG) and C,, (q 


005.00 
O00 

OnelG.919 
0.2498 
Vez 8 
0.3497 
05796 
0.4765 
0.4997 


Z 


woes i 
ys) Sees: 
ooo 70 
oe 2 SZ 
528) SS 1a3) 
woo Ws 
wo oS 
Oo 
oe 7 


Je bOR SIGE S 


C..(q*) 

-00000 
00002 
-00006 
-00010 
-00014 
-00018 
~0G023 
-00031 
-00034 


Slope of Gr. versus q? using Gross' correction is 
_ 00220) Ei). vee 


eope of Gp versus q? UstnG@mearian Ss correction is 
: 0.0214 + 0.0022. 


TABLE D-4 


COMPARISON OF C,.(q*) AND cL (GG ) FOR LF15 


ma 
q* eaae) q? c..(a* ) eaten 
C,(q°) 
0.05 .96900 0.0500 .96900 ~.00000 
O.10 .94002 0.1000 .94003 -.00002 
0.20 .88731 0.1999 .88737 -.00006 
0.25 86324 0.2498 {23922 -.00010 
0.30 .84049 0.2998 .84060 =, OO101.2 
0.35 .81895 0.3497 .81909 -~.00018 
0.40 .79851 0.3996 . 79869 -.00023 
0.477 . 76889 0.4765 .76913 -.00030 
0.50 . 76056 0.4997 .76082 -.00034 


Slope of GED versus q* using Gross' correction is 


Omeogss + 020022. 
pepe OL Gr versus q* using Friar's correction is 
i 0.0180 + 0.0022. 
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APPENDIX E 


COMPUTER PROGRAMS 


Six FORTRAN programs were used to perform calculations 
needed for this work. They are TOPRADCR, TOPROTON, TOPDEUTR, 
BERHFORM, TRADFORM and BERRELCR. The first three were 
written by Topping and can be found in Reference [5]. 
BERHFORM and TRADFORM calculate the experimental cross 
sections for each of the four channels considered for 
hydrogen and deuterium respectively. BERHFORM is repro- 
duced on the following pages. BERRELCR was used to 


compare the relativistic corrections of Gross and Friar. 
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